Creation and Manipulation of Anyons in the Kitaev Model 
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We analyze the effect of local spin operators in the Kitaev model on the honeycomb lattice. We 
show, in perturbation around the isolated-dimer limit, that they create Abelian anyons together 
with fermionic excitations which are likely to play a role in experiments. We derive the explicit 
form of the operators creating and moving Abelian anyons without creating fermions and show that 
it involves multi-spin operations. Finally, the important experimental constraints stemming from 
our results are discussed. 
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Anyons are particles which obey quantum statistics dif- 
ferent from bosons and fermions [1]. In condensed matter 
physics, such weird objects are expected to arise, for in- 
stance, in fractional quantum Hall systems [2], p-wave 
superconductors [3], or Josephson junctions arrays [4], 
but have never been directly observed yet. Within the 
last few years, anyons have also been shown to be of 
primer interest for topological quantum computation [5] . 
Following Kitaev's seminal work related to this problem- 
atics [6] , several spin systems with emerging Abelian and 
non-Abelians anyons have been proposed [7, 8], together 
with experimental proposals susceptible to capture their 
fascinating braiding properties [9, 10, 11]. In this per- 
spective, Kitaev's honeycomb model [7] which only in- 
volves two-spin interactions seems easier to realize than 
the toric code model [G] based on four-spin interactions. 

As shown by Kitaev using 4th order perturbation the- 
ory, the honeycomb model can be mapped onto the toric 
code, thus proving the existence of anyons in the model 
7]. However, orders lower than 4 also bring their physi- 
cal contributions and deserve a careful treatment, espe- 
cially when considering correlation functions and, more 
importantly, the action of spin operators. The aim of 
this Letter is to shed light on these issues and to bring 
some insights for experiments. We show that the original 
spin operators applied onto the ground state create both 
anyons and fermions and we give the explicit form of an 
operator creating anyons without fermions. At lowest 
nontrivial order, this operator involves a superposition 
of states involving one-spin and three-spin operations on 
the ground state. To conclude, we discuss the constraints 
raised by our results for experiments in optical lattices. 

The model — Kitaev's honeycomb (or brickwall) model 
is a two-dimensional spin-1/2 system described by the 
Hamiltonian 

# = - £ £ Jovfof, (i) 

a=x,y,z a — links 

where the erf's are the usual Pauli matrices at site i. 
Without loss of generality [7], in the following, we assume 



J a > for all a and J z > J x ,J y . 

For our purpose, it is convenient to map the spin model 
(1) on the honeycomb lattice onto an effective spin hard- 
core boson problem on a square lattice. This mapping is 
achieved via the following rules 

<.=rUH + h) » <o=iri(bl-b i ), (2) 
<.=r? , < = r?(l-26+& i ), 

where b\ (&^) is the creation (annihilation) operator of a 
hardcore boson at site i of the square lattice depicted in 
Fig. 1 and t" the Pauli matrices of the effective spin at 
the same site. Here, we have attached each site of the 
honeycomb lattice (• or o) to the z-dimer it belongs to. 

In the following, we set J z = 1/2 so that the Hamilto- 
nian reads 

N 

H = -— + Q + T + T +2 + T_ 2 , (3) 
where N is the number of z-dimers, 

Q = £&&> (4) 

i 

T = - £ ( J, t$+ ni + J y t? n > + h.c.) , (5) 
i 

T +2 = -J2( J * + J v 4 +n2 ) = T -2> (6) 

i 

with hopping and pair creation operators t and v 

t ni =bt +ni b i r^n 1 , t n2 = -^1^^1,(7) 
v? ni =bt +ni bUl+ ni , ^=ib\ +n2 bU +n2 Tl.($) 

We emphasize that this mapping is exact and simply 
relies on an alternative interpretation of the four pos- 
sible spin states of a dimer (see [12, 13] for details). 
A key feature of this model is that H commutes with 
W p = ala\alala\al = (-l) fc I fc i+ b DW^ r~ t% r*, with 
notations given in Fig. 1, so that it describes a Z2 gauge 
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FIG. 1: (color online). Mapping of the honeycomb (brick- 
wall) lattice (a) onto an effective square lattice (b) with unit 
basis vectors ni and n?.. The numbering of the sites of a 
plaquette p is shown in both cases. 



theory. Further and more interestingly, it is exactly solv- 
able via several fermionization methods [7, 14, 15]. Its 
phase diagram consists in a gaplcss phase for J x + J y > J z 
and a gapped phase otherwise. The latter contains high- 
energy fcrmions and low-energy Abclian anyons which 
are the focus of the present work. 

Perturbation theory — Here, we focus on the isolated- 
dimcr limit (J 2 3> J x , Jy) and, following Kitaev, we per- 
form a perturbativc analysis of the gapped phase. A 
reliable tool to investigate this problem is the continu- 
ous unitary transformations (CUTs) method [1G], used 
in its perturbative version [12, 17, 18]. Within this ap- 
proach, the Hamiltonian (3) is transformed into an ef- 
fective Hamiltonian H c g which conserves the number of 
bosons ([H e g,Q] = 0). The cornerstone of the method 
is that H c ff is unitarily equivalent to H, i. e., there ex- 
ists a unitary transformation U such that H c s = HU . 
This implies that H and H e g have the same spectrum 
but not the same eigenstates. At order 1, the full effec- 
tive Hamiltonian reads H g = —N/2 + Q + To, and leads 
to high-energy fermionic excitations made of a hardcore 
boson and a spin-string [12, 13]. At order 4, the effective 
Hamiltonian in the low-energy (no fermion) subspace is 
the toric code Hamiltonian [71 



H. 



eff 



N 



J 2 y 



J T + J.. •/, J, 



S 



27V 



E^p- ( 9 ) 



A simple way to obtain this result is to study the ac- 
tion of H, in perturbation, on the low-energy eigenstates 
of the unperturbed Hamiltonian (J x = J y = 0) which 
are ferromagnetic configurations of z— dimers. However, 
we emphasize that the low-energy eigenstates of H are 
not solely built from these states, but involve high-energy 
(antiferromagnetic) configurations, even at order 1 [19]. 
This point is the source of the discrepancy between the 
correlation functions computed in [20] and exact results 
[12, 21]. The perturbative expansion of the exact ground 
state spin-spin correlation functions at order 2 read 



Ko<.) 



1 - {Jl + Jy 



(10) 

(11) 



These results clearly show that sticking to ferromagnetic 
dimcr configurations fails to capture the form of the low- 
energy states [19, 20]. 

The Hamiltonian (9), known as the toric code [6], has 
its ground state in the vortex- free sector (W p = +l,Vp) 
and contains Abelian anyons defined as plaquette exci- 
tations (W p = —1). The real challenge is obviously to 
observe these anyons experimentally. In this perspec- 
tive ultracold atoms in optical lattices are good candi- 
dates [9] although devices using polar molecules have also 
been proposed [10]. Recently, several detection protocols 
in optical lattices based on the possibility to perform 
sequences of single-spin operations have been proposed 
[11, 22]. Unfortunately, although such operations create 
and move anyons, they also create fermionic excitations 
which are susceptible to spoil the detection process. Our 
aim is to show that this crucial fact, always neglected in 
previous studies, is of major importance for experiments. 

Single-spin operations — Let us first discuss the action 
of single-spin operations on the ground state and focus 
on r? = erf.. This operator is used, in Refs. [11, 23], 
as the basic tile to create two anyons on left and right 
plaquettes of the site i [see inset in Fig. 2 (left)]. As 
mentionned above the main problem is that it does not 
simply transform the ground state |0) (belonging to the 
vortex- free sector) into the lowest-energy state of the two- 
vortex configuration one seeks to obtain. To make this 
statement quantitative, we consider the spectral weights 



^ = El^?.P2}>n,fe|Tf|0)f 



(12) 



which obey the sum rule J2 n In = 1- Here, \{p}, n, k) de- 
notes the eigenstate of H in a sector given by an anyon 
configuration W p = —1, and n high-energy quasiparti- 
cles with quantum numbers k. Plaquettes p" and p^ 
are neighbors of i and depend on a = y, z (see insets 
of Fig. 2). With these notations, the eigenstates of H g 
are given by W\{p},n,k). To compute I*, one is thus 
led to determine how the observable r? is renormalized 
under the unitary transformation U. Within the CUT 
formalism, this can be achieved efficiently order by order 
in perturbation [13, 24]. At order 1, one gets 



U\?U = r?[l + (j x vi 



- ni + Jy4- n2 + h.c.)] , (13) 



which clearly shows that pairs of quasiparticles are cre- 
ated or annihilated. In addition, one can see that the 
action of r? on the ground state yields a superposition of 
states with the appropriate two-vortex configuration, but 
with different number of quasiparticles. We computed 
the renormalization of r? up to order 6 but we empha- 
size that quasiparticles arise at order 1 which makes this 
phenomenon strongly relevant. The perturbative expan- 
sion of U'rfU allows a straighforward calculation of the 
overlap probability which reads at order 6: I§ = 1 — i| 



3 




FIG. 2: (color online). Behavior of the spectral weights I„ 
(left) and 1% (right) for fermion numbers n — 0, 2, 4, as a 
function of the coupling J = J x = J y for J~ — 1/2. Grey 
plaquettes in the insets p* and show where anyons 

are created under the action of r? (left) and tV (right). 



with 



yJ x + J y 



- (J 4 

2 1 x 



9 (JxJy 



^J x Jy 



JxJy) 



(14) 



Note that only even n contributions arc nonvanishing. In 
addition, we found that I^>4 = at order 6, showing the 
sum rule is fulfilled. 

Multiple-spin operations — The same approach can be 
used to analyze the action of t? — o\ .erf on the ground 
state, which creates two anyons on the up and down pla- 
quettes of the site i [see inset in Fig. 2 (right)] as well as 
quasiparticles. As previously, we computed the spectral 
weights up to order 6 and obtained: Iq = 1 
with 
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Jijl 



(15) 
(16) 



The behavior of I* and 1% are plotted in Fig. 2 as a 
function of the parameter J x = J y = J. One can di- 
rectly observe the growing weight (for increasing J) of 
the excited states after the operations t z or r y have been 
applied on the ground state. 

Other important quantities are string operators which 
are crucial in experiments to perform braiding of anyons. 
To investigate the effect of such multiple-spin operations, 
let us consider a string 5" = Yl a=1 m of . along a horizon- 
tal line of the original brickwall-lattice (see Fig. 3 with 
m = 3 for notations). When such a sequence is applied 
on the ground state |0), the probability to find the final 
state in the lowest-energy state with anyons at plaquettes 



FIG. 3: (color online). Action, in the vortex-free sector, of 
the string operator S = cr^ .af^.af^.. Each operator flips 
the two plaquettes adjacent to the z-dimer it is attached to 
but also creates fermionic excitations (not shown). 



1 and (m + 1) is given, at order 2, by 

| ({1, m + 1}, 1 5*10) | 2 = 1 — rn( J 2 + J?) 



(17) 



which obviously coincides with /g for m = 1. Note that 
we dropped the k index which is useless for states without 
fermions. This result shows that, at this order, the role 
played by the excited states is directly proportional to 
the length of the string. 

Anyons without fermions — At this stage, we have 
shown that the action of the original spins a" on the 
ground state docs generate unwanted excited states as 
can be seen at lowest nontrivial order in the perturbation 
theory. One may wonder how to create anyons without 
fermions. In other words, we wish to determine the op- 
erators that only flip two W p 's and nothing else. 

As an example, we focus on the operator creating two 
vortices at the left and right plaquettes of a given site i 
[see inset Fig. 2 (left)]. Let us denote this operator fij, 
and show how to compute its perturbative expansion, 
Hi = X^fceN ' wnere ^i fc ' contains all operators of or- 
der k and thus associated to J l x J™ (with l + m = k). The 
aim is to obtain £li e ff = l/iQiU = r? which indeed leads 
to Iq = 1. At order 0, operators are not renormalized so 



(0) 



The renormalization 



that one obviously has Q.\ 
of fij under the unitary transformation U reads 



(fc),[Z] 
*eff 



(18) 



feeN 



ken let 



where Cli^f is of order (k + I). Order being already 
fulfilled, let us write this constraint at order 1, using 



Eq. (13) and ft;™ 



(A') 



o = a 



(0), 



(1).[0] 



(Jxvt- ni +Jyvt- n ,+ h.c.)+0« 

(19) 



This equation yields so that using the inverse map- 
ping of Eq. (2) one finally gets, in the original spin lan- 
guage and at order 1 



ft* 



■A 



(20) 



Jx \ (J i-n 1 ^ (7 i-Ti 1 ,o (7 i^ 



J V [ <J i-n 2 ^ a i-n 2 ,o <J i,' 



v y z 
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This expression shows that to create anyons without 
fermions, one must be able to build a delicate superpo- 
sition of single spin- flip and three spin- flip states. Fur- 
thermore, the weights of the latter have to be fine-tuned 
since they depend on the precise values of the couplings 
which may be experimentally challenging. Note also that 
higher order corrections would even involve more spins. 

Experimental discussion — Let us now put our results 
in an experimental perspective, setting, for simplicity, 
J x = J y = J and 7 = J/J z . Equation (17) shows 
that the repeated action of m spin-flips decreases the 
weight of the (pure) two-anyon state in the low-energy 
subspace. One should thus work with small enough val- 
ues of 7. For a reasonable number m ~ 25 of operations 
needed to perform a braiding of anyons, and assuming 
that m7 2 /2 ~ 1/2 [25] leads to conclusive experiments, 
one is led to choose 7 < 0.2. 

From a practical point of view, we also wish to mention 
that detecting anyons using ground-state two-spin corre- 
lation functions after an anyonic braiding is rather diffi- 
cult. Indeed, contrary to what is mentioned in [11], these 
functions do not change sign in the presence of anyons. 
The latter are only responsible for 4th order corrections 
in (erf cf.) whose leading term is of order [12], and 
for 3rd order changes in « i0 < +nii .) and (<7? aj' +n2i .) 
whose leading terms are of order 1 [13]. In this respect, 
the more efficient and sophisticated set-up recently pro- 
posed in [22] is an interesting alternative. 

We end this section with some remarks about tem- 
perature issues in experiments, based on the orders of 
magnitude J z /h = 5kHz and J/h = 1kHz [26] (7 = 0.2). 
To perform anyonic interferometry experiments [11, 22], 
the temperature of the system should be small enough 
to prevent thermal excitations of unwanted anyons. As- 
suming thermal equilibrium, the typical low-energy scale 
J off = J z 7 4 /16 appearing in front of the plaquette term 
in Eq. (9) leads to the constraint T <C J c ffA*B — 20pK. 
Such a temperature is lower than recently reached tem- 
peratures [27], but this problem may be circumvented by 
working out of equilibrium. 

To sum up, we have shown that unless experimentalists 
manage to apply the operators given in Eq. (20), they will 
have to make sure the constraint 7 = J/J z < 0.2 is ful- 
filled. This constraint, arising from the renormalization 
of spin operators at order 1 in perturbation, could lead to 
a problematic upper bound on the temperature. The lat- 
ter indeed scales as 7 4 since the toric code Hamiltonian 
arises as an effective Hamiltonian of Kitaev's honeycomb 
model only at order 4 in perturbation. Finally, let us 
mention that the time evolution of the system between 
spin operations may also lead to an unwanted dynamics 



of the fermions but this difficulty can be circumvented by 
freezing the system via an increase of the optical lattice 
barriers. 
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